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Abstract 

We investigate observational constraints on the variable Chaplygin gas model from the 
gold sample of type la supernova data and the recent measurements of the X-ray gas mass 
fractions in galaxy clusters. Combining these databases, we obtain a tight constraint on 
the two model parameters. Our results indicate that the original Chaplygin gas model is 
ruled out by the data at 99.7% confidence level. 
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1 Introduction 



Recent observations of type la supernovae (SNe la) suggest that the expansion of the 
universe is accelerating and that two-thirds of the total energy density exists in a dark 
energy component with negative pressure £Q. In addition, measurements of the cosmic 
microwave background (CMB) j2] and the galaxy power spectrum [2j also indicate the 
existence of the dark energy. The simplest candidate for the dark energy is a cosmological 
constant A, which has pressure P\ = —p\. Specifically, a reliable model should explain why 
the present amount of the dark energy is so small compared with the fundamental scale 
(fine-tuning problem) and why it is comparable with the critical density today (coincidence 
problem). The cosmological constant suffers from these problems. One possible approach 
to construct a viable model for dark energy is to associate it with a slowly evolving and 
spatially homogeneous scalar field 0, called "quintessence" |U |3] . Such a model for a broad 
class of potentials can give the energy density converging to its present value for a wide set 
of initial conditions in the past and possess tracker behavior. 

Recently, the Chaplygin gas model was proposed as an alternative to the cosmological 
constant in explaining the accelerating universe jH] . The Chaplygin gas is characterized by 
an exotic equation of state P = —A/p, where A is a positive constant. An attractive feature 
of the model is that it can naturally explain both dark energy and dark matter. The reason 
is that the Chaplygin gas behaves as dust-like matter at early stage and as a cosmological 
constant at later stage. Some possible motivations for this model from the field theory point 
of view are discussed in Refs. [7j. The Chaplygin gas appears as an effective fluid associated 
with (i-branes [H] and can also be derived from the Born-Infeld action [Oj. An interesting 
range of models was found to be consistent with SN la data ^0], CMB experiments JTI] 
and other observational data [TJ]. The Chaplygin gas model has been extensively studied 
in the literature [T3] . 

However, the Chaplygin gas model produces oscillations or exponential blowup of the 
matter power spectrum that are inconsistent with observation ^3]. In Ref. ^H], we con- 
sidered a variable Chaplygin gas (VCG) model and showed that it interpolates between 
a universe dominated by dust and a quiessence-dominated one described by the constant 
equation of state. Furthermore, we showed that the model corresponds to a Born-Infeld 
tachyon action ^Hj. Recently, the model parameters were constrained using the location of 
peaks of the CMB spectrum and SN la data (TH]. In this paper we consider observational 
constraints on the VCG model from the gold sample of 157 SN la data and the recent mea- 
surements of the X-ray gas mass fractions in 26 galaxy clusters. We perform a combined 
analysis of these databases and obtain the confidence region on the two parameters. Our 
results indicate that the original Chaplygin gas is ruled out at 3a confidence level. 
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2 Variable C haply gin Gas 



Let us now consider the VCG characterized by the equation of state 

Pv 

where A(a) is a positive function of the cosmo logical scale factor a. This assumption is 
reasonable since A(a) is related to the scalar potential if we take the Chaplygin gas as a 
Born-Infeld tachyon field [01 • In a spatially flat Friedmann- Robertson- Walker (FRW) 
universe, the energy conservation equation is 

dpv g Pv + Pv 

da a 

By inserting Eq. into the above equation, one finds that the VCG density evolves as 

1/2 



p v (a) = a 3 



6 J A(a)a 5 da + B 



(3) 



where B is an integration constant. Given a function A(a), Eq. © allows us to obtain a 
solution p v (a). We assume A(a) is of the form A(a) = A a~ n , where A Q and n are constants. 
This ansatz has the following important features: a) the VCG model with n = reduces 
to the original scenario, b) we will see that the VCG behaves as a quiessence rather than a 
cosmological constant at late times and c) the function p(a) can be calculated analytically. 
Then from Eq. J3J) it follows that 



6 A B . lS 



n a' 1 ay 



Note that n = corresponds to the original Chaplygin gas scenario, in which the Chaplygin 
gas behaves initially as dust-like matter and later as a cosmological constant [Hj. However, 
Eq. (jlj) shows that, in the VCG scenario, it interpolates between a dust-dominated phase 
and a quiessence-dominated phase described by the constant equation of state w = — 1 + 
n/6 jUj. From Eq. (J3| we get the present value of the energy density of the VCG 




Prt = \* ^o + 5, (5) 



where the present value of the scale factor is normalized to unity, i.e., ao = 1. Defining 
B s = B/pl , Eq. (|4|) takes the form 

p v {z) = p v0 [B s (l + zf + (1 - B s )(l + z) n ] 1 ' 2 , (6) 

where z = 1/a — lis redshift. In the spatially flat FRW metric the Friedmann equation 
can be written as 



H 2 = j( Pb + Pv ), (7) 



2 



2 



where H = a/a is the Hubble parameter, k 2 = 8ttG is the gravitational coupling and pb is 
the energy density of the baronic matter. Substituting Eq. (jHJ) into the Friedmann equation 
(J7J) gives 



H\z)/Hl = Q b (l + zf + (l-Q b ) 

x 'B s (l + z) 6 + (l-B s )(l+z) 
= E 2 (zMh, B s ,n), 
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where H = 100/ikms _1 Mpc _1 is the present value of the Hubble parameter and is the 
density parameter of the baryonic matter component. Then it is straightforward to show 
that the luminosity distance di and the angular diameter distance d& in the spatially flat 
FRW universe are respectively given by 



d L 



:i+z) 



Ho 

(1 + z)~ 2 d L . 



dz 
o Wz) 



(9) 
(10) 



3 Supernova la Constraints 

Let us now consider constraints on the VCG model from the gold sample of 157 SN la data 
compiled in Ref. ^H]- The parameters in the model are determined by minimizing 

2 _ [/^mod(-^i) Oft, B s , Tl) /^obs(^i)] ('\'\\ 

1=1 

where Oi is the total uncertainty in the observation, /x b s is the observed distance modulus 
of SNe la, and fi m od( z i) is the theoretical distance modulus 

Mmod(^)=51og 10 ^^ + 25. (12) 
Mpc 

To determine the likelihood of the parameters B s and n, we marginalize the likelihood 
function L = exp(— % 2 /2) over h and We adopt Gaussian priors such that h = 0.72±0.08 
from the Hubble Space Telescope Key Project ^Hl and VL^h 2 = 0.0214 ± 0.0020 from the 
observed abundances of light elements together with primordial nucleosynthesis j20] • The 
results of our analysis for the VCG model are displayed in Fig. [T] We show 68.3%, 95.4% 
and 99.7% confidence level contours in the (B s , n) plane. The best-fit model parameters 
and marginalized la error bars are B s = 0.223lo;o59 an d n = —3-0-l'.i with xLin = 174.253. 
The results show that the two parameter B s and n are highly degenerate. It is interesting to 
note that the dark energy component with w < — 1 is favored, which allows the possibility 
that the dark energy density in increasing with time. 
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Figure 1: Probability contours for B s versus n are shown at la, 2a, and 3a when f2 tot = 1. 
These constraints use the gold sample of 157 SN la data 
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Figure 2: Probability contours for B s versus n are shown at la, 2a, and 3a when f2 tot = 1. 
These constraints use the X-ray gas mass fraction in 26 galaxy clusters 
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4 Constraints from the X-ray Gas Mass Fraction of 
Galaxy Clusters 

The matter content of rich clusters of galaxies is thought to provide a fair sample of the 
matter content of the universe as a whole. The observed ratio of the baryonic to total mass 
in clusters should therefore closely match the ratio of the cosmo logical parameter Q b /Q m , 
where fit, and fi m are the mean baryon and total mass densities of the universe in units of 
the critical density. The combination of robust measurements of the baryonic mass fraction 
in clusters with accurate determinations of f2& from cosmic nucleosynthesis calculations 
can therefore be used to determine Q m [2I|. The measurements of the apparent redshift 
dependence of the baryonic mass fraction can also, in principle, be used to constrain the 
geometry and dark energy density of the universe |22j- The first successful application 
of such a test was carried out by Allen et al. using a small sample of X-ray luminous, 
dynamically relaxed clusters with precise mass measurements, spanning the redshift range 
0.1 < z < 0.5 [2H] (see also |24j). Recently Allen et al. present a significant extension 
and obtained a tight constraint on the mean matter density and dark energy equation of 
state parameter. The clusters sample is significantly larger and includes 26 X-ray luminous, 
dynamically relaxed systems spanning the redshift range 0.07 < z < 0.9 (22] • We will use 
this database to constraint the VCG model. To determine the confidence region of the 
model parameters, we use the following \ 2 function 



26 \f mod ( z -h fL B n) — f obs 

2 V"^ L gaS •'gas.*] /-|o\ 

Xxray = ~2 ' 



where /^f f is the measured X-ray gas mass fraction f gas with the defaut standard cold 
dark matter (SCDM) cosmology, o" gaSj j is the symmetric root mean square errors, and 
is the model function 



J gas W (1 + 0.19^/2)^ 



h df™( Zl ) 

0.5 dT d {zi) 



3/2 



(14) 



where the bias factor b = 0.824 ± 0.089 25. is a parameter motivated by gas dynamical 
simulations, which suggest that the baryon fraction in clusters is slightly depressed with 
respect to the universe as a whole, and the effective matter density parameter is 

a£ = n b + (i-n b ),JF a . (15) 

Adopting Gaussian priors such that h = 0.72 ± 0.08 and Vt b h 2 = 0.0214 ± 0.0020, the 
68.3%, 95.4% and 99.7% confidence level contours in the (B s , n) plane are shown in Fig. EJ 
The best-fit model parameters and marginalized la error bars are B s = 0.049^015 an d 
n = 0.5lj;;§ with xLm = 24.437. The results favor the original Chaplygin gas model with 
n = 0. 
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Figure 3: Probability contours for B s versus n are shown at la, 2a, and 3a when f2 tot = 1. 
These constraints use the gold sample of 157 SN la data JHI and the X-ray gas mass fraction 
in 26 galaxy clusters |2H] . 



5 Combined Analysis, Conclusions and Discussions 

Let us now consider a combined analysis of the constraints from the gold sample of SN la 
data and the measurements of the X-ray gas mass fraction in galaxy clusters. The fit is 
done by minimizing the following \ 2 function 

Xtot = XSN + xlray ■ (16) 

Fig. El shows the constraints on B s and n obtained from the analysis of the combined the 
SN la and the X-ray gas mass fraction data set. We see that SNIa+/ gas data set provides a 
remarkably tight constraint in the (B s , n) plane, with best fit values B s = 0.055i !oi2 (1°" 
error bar) and n = 1.70lto[ii (lc error bar) with Xmm = 205.457. We find that the original 
Chaplygin gas model is ruled out by the data at 99.7% confidence level. Fig. |3]and Fig. |5] 
show the probability distribution for n and B s , respectively, marginalized over B s and n in 
the spatially flat FRW universe. 

The VCG model, which plays the role of both dark energy and dark matter in terms 
of a single component, deserves to explore its various observational effects ^Hj. In this 
paper we have focused our attention on the gold sample of SN la data and the recent 
measurements of the X-ray gas mass fractions in galaxy clusters. Adopting simple Gaussian 
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Figure 4: Probability distribution for n marginalized over B s obtained from the analysis of 
the combined SNIa+/ gas data set in the VCG model. 




Figure 5: Probability distribution for B s marginalized over n obtained from the analysis of 
the combined SNIa+/ gas data set in the VCG model. 
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priors of h and Qbh 2 , we have obtained stringent constraints on the two parameters, B s 
and n, which describe the effective density parameter of the dark matter component and 
the effective equation of state of the dark energy component, respectively. The combined 
analysis of these databases shows that the original Chaplygin gas model is ruled out and the 
universe tends to be quiessence-dominated rather than phantom-dominated at 3a confidence 
level in the VCG scenario. It would be interesting to investigate the evolution of density 
perturbations in this VCG model. 
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